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Abstract 

The exact renormalization group equation for pure quantum gravity is derived for an 
arbitrary gauge parameter in the space-time dimension d = 4. This equation is given by a 
non-linear functional differential equation for the effective average action. An action func- 
tional of the effective average action is approximated by the same functional space of the 
Einstein- Hilbert action. From this approximation, /3-functions for the dimensionless Newton 
constant and cosmological constant are derived non-perturbatively. These are used for an 
analysis of the phase structure and the ultraviolet non-Gaussian fixed point of the dimen- 
sionless Newton constant. This fixed point strongly depends on the gauge parameter and 
the cutoff function. However, this fixed point exists without these ambiguities, except for 
some gauges. Hence, it is possible that pure quantum gravity in d = 4 is an asymptotically 
safe theory and non-perturbatively renormalizable. 
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1 Introduction 



As is well known, quantum gravity (QG) must be treated non-perturbatively. This is because, 
in the space-time dimension d = 4, the L-loop perturbative calculations in Einstein gravity 
cause divergences that are proportional to the L + 1 powers of the curvature tensor. Hence, the 
renormalization of these terms requires infinite number of couplings. Thus, QG is called (per- 
turbatively) non-renormalizable. However, if QG is an asymptotically safe theory, it becomes 
(non-perturbatively) renormalizable 0]. An asymptotically safe theory is specified by the exis- 
tence of a ultraviolet (UV) non-Gaussian fixed point (NGFP). The asymptotically safe nature 
of QG is suggested by d = 2 + e gravity theory It is also expected that this nature will be 

maintained in d = 4. However, an ordinary perturbative e-expansion is an asymptotic expan- 
sion. Thus, the large order behavior of e is not reliable. Hence, to guarantee the existence of the 
UV NGFP in more higher dimensions, the appropriate method such as the Borel resummation 
must be applied. 

The exact renormalization group equation (ERGE) || used in this article is one of the non- 
perturbative methods in field theories. In addition, the applicability of the ERGE exceeds the 
e-expansion and the 1/iV-expansion In the previous article []|], we used the formulation of 
the ERGE for pure QG f§, and clarified that QG has the UV NGFP in 2 < d < 4. This result 
suggests that QG is an asymptotically safe theory and (non-perturbatively) renormalizable. 
However, this UV NGFP depends on the gauge parameter and the cutoff function. Hence, the 
purpose of this article is to study the effect of that dependence to the UV NGFP. 

The cutoff function dependence is an artifact problem. The origin of this problem is stated 
as follows. The ERGE is formulated as a non-linear functional differential equation for the 
effective average action. The effective average action is an infrared (IR) cutoffed 1PI effective 
action, and corresponds to a coarse grained free energy. To formulate the effective average 
action, we must introduce the cutoff function. As we will discuss in Sec. 2, the profile of the 
cutoff function is arbitrary. If we can solve the ERGE without any approximations, this problem 
will not appear m. However, it is impossible to solve the ERGE without any approximations. 
Hence, to reduce this equation to the calculable form, we must truncate the functional space of 
the effective average action. For this truncation, the non-linear functional differential equation 
is reduced to a set of the non- linear differential equations. This truncation causes the cutoff 
function dependence. 

The origin of the gauge dependence in the ERGE is same as that of the ordinary field theory. 
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If we only discuss the existence of the NGFP, this problem is not so serious. This is because, 
the non-physical quantity such as the /3-function and the fixed point (FP) may depend on the 
gauge. However, if this dependence causes the disappearance of the NGFP, it becomes problem, 
because the gauge dependence changes the phase structure of the theory. We have been believed 
that Abelian and Yang-Mills gauge theories have only the GFP. Thus, the gauge dependence of 
the NGFP have not been discussed. However, the existence of the NGFP is an important point 
in QG. Though we expect that the NGFP should remain for any gauges, a strict proof is not 
known. 

This paper is organized as follows. In the next section, the formulation of the ERGE for 
pure QG is reviewed (for details see In Sec. 3, the functional space is approximated by 
the same space of the Einstein-Hilbert action. In there, the /3-functions for the dimensionless 
Newton constant and cosmological constant including the constant gauge parameter are derived 
by a slightly different method from other formulations || ||, |l0|. In Sec. 4, these /3-functions are 
used to study the effect of the gauge and cutoff function dependence to the UV NGFP. Section 
5 is devoted to summary and discussion. 

2 Formulation and approximation of the ERGE 

To derive the ERGE for pure QG, we define the scale dependent generating functional W k of 
the connected Green functions. Here, k is the IR cutoff scale in the Euclidean momentum space. 
Hence, W k is the IR cutoffed generating functional and defined by 

e W k [t,a,a;/3,r;g) = J J)hVCVC exp {- S grav [j] - S g . L [h; g] 

-S F .p. [h, C, C; g] - S e . s . [t, a, a; (3, r; g] 

-A k S^ v [h;g}- A k S^[C,C;g]} . (1) 

Here, SgravM is a general functional of a quantum metric 7^. In the background field method, 
the quantum metric is decomposed as the background metric g^ and the fluctuation from the 
background h^ v . The functional S g ^\h;g] in Eq. ([!]) is the gauge fixing term given by 

where a is a gauge parameter. In the harmonic gauge = 0, the explicit form of F^ is given 

by 
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where « = (32vrG)~ 1/2 and G is the bare Newton constant. The covariant derivative D is 
constructed from g^ v . In Eq. ([[]), Sp,p\h,C,C;g] is the Faddeev- Popov ghost term given by 

S F . P .[h,C,C;g] = ~V2j dfx^C^M^gf V C V , 

where C M and C p are the ghost and anti-ghost fields respectively. Here, 

M[ T ,g}\ = r p 9 aX Dxh P ,D a + lav D p ) - g pa r X Dx lau D p . 
The functional S e , a , in Eq. (|l|) is the external source term and given by 

Se. s . [t, a, a; p,r;g] = -J d A x^g {iTh^ + + a^C p 

+P flu (s B h flu ) + T fl (s B cn}- 

Here the sources (3^ v and r p couple to the BRST variations of h^ v and respectively. 

In Eq. (|l]), AfcS' grav [/i; g] and AfcS' gh [C, C; g] are the cutoff actions. These terms control the 
propagation of fields to have the momentum k < p < ko. Here ko is a UV cutoff scale of the 
theory. Hence these are given by a quadratic form of and the ghost fields. The explicit 
forms of these terms are given by 

A k S^[h;g] = \k 2 f d d x^gh, v {Rr v [9]r pu K^ (2) 
A k S^ h [C,C;g] = V2 J d d x^C p Rf \g]C, (3) 

where the cutoff operator -Rf rav [g] and -Rf h [<?] are defined by 

Rf av [g] = {Zf™Y UP<J k 2 R {Q) {-D 2 /k 2 ), (4) 
Rt[9] = Zl h k 2 R(°\-D 2 /k 2 ). (5) 

Here, (Zf * r )i»'P° and Zf is the renormalization factor of h pv and the ghost field respectively. 
A convenient choice of the cutoff function R^\—D 2 /k 2 ) is 

1 + e -D 2 /k 2 _ e -D 2 /k% 



R W(-D 2 /k 2 ) 



e I ~o — e 

If we take the limit ko — > oo, we have 



R M(-D 2 /k 2 ) 



-D 2 /k 2 



e -D 2 /k 2 _ i 

Hence, the constraints of the cutoff function R(°>(u) are given by 



lim R.W(u) = 1, lim R^(u) = 0. (6) 

u^O «^oo 



Any functions satisfying Eq. (q) are applicable. Hence, the profile of the cutoff function is 
arbitrary. In subsection 4.2, the effect of this ambiguity to the UV NGFP will be considered. 
The effective average action is defined by 



rfc[/i,£,£;/?,r;#] 



d*Xy/§ {t^ V + a„e + ^} ~ W k [t, a, a; (3, r; g] 
-A k S^[h-,g]-A k S*[S,S;g]. 



(7) 



Here, scale dependent classical fields are given by 

1 5W k ^ 1 SW k 



h 



i sw k 



In addition, the classical field corresponding to 7^ is introduced as 

The effective average action has two boundary conditions. One is given by 

lim Tk[h,C,C;/3,r;g] = T 1P1 [h,£,£; f3,r;g], 

in the limit k — ► 0. This is because, all quantum corrections are included in T k= $ in this limit. 
Thus this is equivalent to the ordinary 1PI effective action. The other is given by 



lim T k \g,t,£;P,T;g] 

k—>oo 



Sgravb] + S g s.[g - g;g] + 5 F .p. [g - g,£,£;g] 



(8) 



in the limit k — > 00, Here, we denote huu as g^ v — g^ v . Equation (||) means that T k is coincide 
with the bare action in this limit, since there are no quantum corrections. 

If Eq. (||) is differentiated with respect to t = Ink, and Legendre transformed by Eq. ([/]), we 
obtain 

1. 



d t T k 



-Tr 



(if } + ^R g r)- hh (d t ^Rr v r upa 



1. 



-Tr 



r k 2) +V2Rt - (rf )+V5flf yj] (<W5flf 



(9) 



(2) 

This is the ERGE for pure QG. Here T k is the Hessian of T k with respect to the subscript. 

Though Eq. (|9|) is non-perturbatively exact, the manifest BRST invariance is broken by the 
explicit IR cutoff. To see this case more detail, we consider the Ward-Takahashi (WT) identity. 
Now, the WT identity is given by 

= (S B S M . + S B A k S^ av + d B A k S^ h 



If this is written in terms of the effective average action, we have 

/ A 7s {StM + §3 = » (*™ • <10) 

where F ; k = T k — 5 g .f. is introduced. In the usual field theories, the RHS of Eq. (|l0|) equals 
to zero. In the present case, the existence of the cutoff action makes it proportional to cutoff 
operators. The RHS of Eq. (|l0|) goes to zero in the limit k — > 0, because cutoff operators goes 
to zero in this limit. Hence, an usual field theory is recovered. However, Y k does not disappear 
in the intermediate scale k. Thus, the BRST symmetry is broken in this scale. 

Now, to get the BRST invariant RG flows, we approximate the functional space. As a first 
step approximation, we neglect the evolution of the ghost action and the external source fields. 
From this approximation, the effective averaged action is expected as 

Fk[g,€,€;(3,T;g] = T k [g] + t k [g,g] + S g _t[g - g;g] + S F .p.[g - <?,£,£;<?] 

-J ^vsir^y+r^iB^)}. (ii) 

Here Sg.f. and Sf.p. have the same form as in the bare action. The coupling to the BRST 
variations also has the same form as in the bare action. The remaining term is decomposed into 
T k [g] and T k \g,g]. Here, T k [g,g] contains the deviations for g^ v ^ g^ u , and satisfies Tk[g,g = 
g] = 0. The approximated effective average action given by Eq. ( pd] ) satisfies the boundary 
condition of Eq. (||), if these terms satisfy 

lim r k [g] = Sgravb], lim t k [g,g] = 0. 

fe— too fe^oo 

These conditions suggest that setting T k [g,g] = for all k is the candidate to get the BRST 
invariant RG flows. Substituting Eq. Q) into Eq. ([Tq| ), we obtain 

/ ddxC ^ 6 ^M = ~ Yk { R r[9lRf[g]) , (12) 

where means the Lie derivative with respect to The RHS of Eq. (^) is regarded as the 
higher loop corrections if it is evaluated perturbatively (for details see Q). Hence to neglect Y k 
is acceptable in the first approximation. This is consistent with setting T k = for all scales. 
These approximation means that the RG flows are projected on g^ = g^ v in all scales. In the 
background spaces the BRST invariance is preserved. Hence, the projected RG flows moving in 
the background space are regarded as the BRST invariant. 



If Eq. (11) is inserted into Eq. (|9|), the approximated ERGE becomes 

d t T k [g;g]=Sr-St, (13) 
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where, 

T k [g;g]=T k [g]+t[g;g] + S g . { .[g-g;g]. (14) 
This has the boundary condition 

lim T k [g; g] = S ffSV \g] + S g , L [g -g;g}. 

K—>OQ 

In Eq. (Jl3|) , 5| rav and 5| h correspond to the gravitational sector and the ghost sector respectively, 
and are given by 



^grav 



-Tr 



-Tr 



(V 2 ri 2) [^] + Rrmy 1 (d^rm) 



-M&gj + RtlgjY'idtRtlg]) 



(15) 
(16) 



In Eq. ([H]), T^\g;g\ represents the Hessian of T k [g;g] with respect to g^ u at fixed g^ v . 



3 Einstein-Hilbert truncation 

In below we consider the case d = 4. Now, to make problems easier, we truncate the functional 
space of T k [g; g\. The most naive truncation is to take the functional space as a same space of 
the Einstein-Hilbert action. Thus the bare action is given by 

1 



Sgrav [<?] 



WttG 



d i x^g-{-R{g) + 2\} 



where A is the bare cosmological constant. Now, we define the scale dependent couplings as 

G — >■ G k = Zjj\G, A —¥ X k , a — > a k = Zjj\a. 

Hence, from Eq. (|I4|), T k [g;g] is expected as 

T k [g;g] = 2 K 2 Z Nk J d i x^{-R(g) + 2\ k } 

+^Z Nk J d^x^gt v ^fg al3 ){T^g p(J ). (17) 

Here Z^ k is a renormalization factor. If Eq. (|l~7| ) is differentiated with respect to t and projected 
on g^y = g^ u , we have 

d t T k [g;g] = 2k 2 [ d 4 x ^ [-R{g)d t Z Nk + 2d t (Z Nk X k )] • (18) 



This is the LHS of Eq. (jT3|) . Though the differentiation with respect to t brings the term that 
is proportional to the ghost action, this term disappears on g pu = g pv because Tff 'g a f3\g=g = 0. 
Thus the gauge parameter is treated as a constant. This is a problem of the present formulation. 
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Next step is to get the RHS of Eq. @. Firstly, we calculate St given by Eq. (H). Now, 
we can naively write Eq. ( |T5| ) as 



^grav 



2 



(V^Wl +RrMy 1 (K- 2 d t T®[g,g 



(19) 



The FP action T*[g; g] satisfies dtT*[g; g] = 0. Hence, if we are interested in only the FP solution, 
we can neglect the second term in the RHS of Eq. (|l9|). Thus, we have 



^grav 



-cUn/f 



grav r 



where, 



4 grav [ 



Vh, u exp {-Tr d [h; g) - A k S^ v [h; g]} . 

Here, A k S gTa,v is given similarly to Eq. (|2|) except for the change of the field: h 
r^ uad [/i,^] is defined by 

K- 2 T k [g = g + h;g]= T k [g] + 0(h) + Tf^h-g] + 0(P). 

The explicit form is given by 



(20) 
(21) 



h^, and 



Tl uad [h;g] = Z Nk I d 4 x^/g h 



l pV 



per 1 per 



h p(T 



1 - Z Nk J d^g gT (FfKp) {^hpa) , (22) 



where, 



V p 

a p 



u^ pa = K^ p AR-2\ k ) + \(r v R P . + g P aRn-\{R v p^ + R 
-\ (5»R\ + 8»R% + 5 v p R^ + %R" p ) . 

Hence, if the one-loop effective action /f rav is calculated, we can get 5f rav from Eq. (p0|). Now, 
to calculate the one-loop effective action, we decompose h^ iu as 



h 



{Dpil 



1 



pr \^p^u + D vip ) + ( D p D u ~ -^puD j rr ■ -g,,„0. 



1 



(23) 



where (j) = g^ v h^ v , and £j satisfies D^^ p = |^, pT[| . In addition, h pu satisfies D^h pi> = and 



g» v h 



pv 



0. The resulting measure and jacobian are given by 



Vh 



pv 

J 



Vh^V^Va[detJ] 



1/2 



:R 



Ar 



R ® A (0) 



(24) 
(25) 



Furthermore, we take the background as the maximally symmetric. In this background, the 
Riemann tensor R^ up(T and the Ricchi tensor are given by 

1_ - 

Rpu — ~^)ixvR- 

Here, the scalar curvature R characterizes the space. Now, we introduce the constrained opera- 
tors as 

Ao(X)0 = {-D 2 + X 
A lfll/ (X)& = (-D 

From these manipulations, Eq. becomes 



2 ^ u + 6£S£X) h 



^quad 



[kg] = z Nk J d*x^^h^ v A 2 U2\ k + ^h x ^ 



a p 
3(q - 3) 
16a 
, 2{a-l) & 



-2a\ k + ^R)t» 



4a 
a — 3 



a — 1 

a — 3 



+ 



a — 3 
3a — 1 



a^Ao (0)^/Ao (-R/3)<j> 
2a 



3(a - 3) 



3a — 1 



Here, we introduced 



Ai 



R)&, a = w Aq(0)w A, 



1 



(26) 



(27) 



The changes of variables in Eq. (|27| ) cancel the jacobian in Eqs. (24) and (|25|). 

To get the one-loop effective action given by Eq. (pl|), we decompose /i M „ in AjfcS' grav [/i; 5]. Up 
to now, the form of -2?f rav in Eq. (|j) is not specified. Now, we consider i?| rav as the non-trivial 
projection operator which makes A k S srav [h; g] to 



A k S^[h,g] 



Z Nk / (fx^Jg 



+ 



3(a - 3) 

16a 
2(a - 1) 
a — 3 



A (0WAo(-.R/3) 



(28) 
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Here, we introduced the shorthand notation R^°\—D 2 /k 2 ) = §. Hence, from Eqs. ( |26| ) 
and (|2^), we can calculate I^ av by the Gaussian integral. However, that result includes the 
additional zero- modes introduced by the decomposition in Eq. (p3|). Then to remove these 
modes ||, 11 1, we introduce unconstrained operators as 

detAs(X) = detAo(X), 



det A v (X) 

det A T {X) 
Hence, we have 



det Ai (X) det A ( X - -R ) , 



det A 2 (X) det A a ( X - —R ) det A ( A' 



-R 



detZ Nk A T ( k 2 R {0) - 2\ k + -R 



det Z Nk A y ( k 2 R^ - 2aX k + 



2a - 1 



R 



det Z Nk A v ( k 2 R^ - 2\ k + ^R 



detZjvfcAg (k 2 R^ - 2A fc 



Therefore, from Eq. (|2l| 



^grav 



1 



Tr T 



-It, 



.V f A + 



+ -Tr v 
+ ^Tr s 



AT [A a + ^—^R 



MA 



-1 



where, .4, A a and are given by 

A = - D 2 + e R (v)(- D 2/ k 2)- 2 \ k , 
A a = -D 2 + k 2 R (0) (-D 2 /k 2 ) -2a\ k , 



Af 



1 _ 1 ) k 2 R^(-D 2 /k 2 ) + D 2 R^'(-D 2 /k 2 ). 



(29) 



Here, the prime means the differentiation with respect to the argument. The anomalous dimen- 
sion r\ is defined by rj = —dt \n.Z^ k . In Eq. (|29"|), = is taken. In below, we omit the bars 
from the metric and the scalar curvature. 

The remaining term of the RHS of Eq. (|l3|) is 5| h . For the Faddeev-Popov ghost term, we 
do not take into account the renormalization of these field. Hence, i?| h = k 2 R^\—D 2 /k 2 ). In 
the present background, AA = —D 2 — R/A. Therefore, 



-TV V 



M, ( A - -R 



-1 



(30) 
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Here, Aq and Mo in Eq. (|3^) are denned similarly to A and M except for X k = and r\ = 0. In 
below, we denote the RHS of Eq. ( |l3| ) as 5r = 5| rav + 5| h . 

Now to get the coefficients of yfg and y/gR, we expand 5r in terms of the scalar curvature 

R, 



Sr 



^Tr T 
+ ^Tr s 



MA' 1 
MA' 1 



MA- 1 



MA- 1 



Tr 



f,'^Tr T 



v 



+ 



2a - 1 



AM; 



-T^Try 



MA- 



Try 



AAqV \ + 0(R 2 ). 



(31) 



Furthermore, to calculate the traces in Eq. (|3l|), we use the heat kernel expansion: 



W(-D 2 ) 



(47r)- 2 t ri (7)|Q 2 [W] J d A x^g- 
+\Qx[W] J d 4 x^R + 0(R 2 )Y 



(32) 



where, I is a unit matrix and j = T, V, S mean the tensor, vector and scalar respectively. In 
Eq. (|32|), trj(J) simply counts the number of independent degrees of freedom of these quantity: 

tr T (/)=9, tr v (/)=4, tr s (J) = 1. 

In Eq. ([?!), Qt[W], [i = 1, 2) is the Mellin transform of W, and given by 



Qo[w] 

Qn[W] 



W(0), 

1 



dzz n ~ 1 W{z), (n > 0). 



9t (ZNk^k) 



dtZNk 



T(n) Jo 

Therefore, if we insert Eq. ( |32[ ) into Eq. (|3~l|), and compare it with Eq. (|l8|), 

4^(£) [ 6$ 2(-2A fc A 2 ) +^ 1 {-2a\ k /k 2 ) 
-r, {M>\(-2\ k /k 2 ) + 2¥ 2 (-2a\ k /k 2 )}] , 

-30<5>l(-2\ k /k 2 ) - 6(2q - l)$ 2 2 {-2a\ k / k 2 ) 
-^{{-2\ k /k 2 ) + 4<S>\(-2a\ k /k 2 ) - 12*|(0) - 8^(0) 
-77{-15$l(-2A fc /fc 2 ) - 3(2a - l)$ 2 (-2aA fc /£; 2 ) 



(33) 



1 ^ 
12k 2 (4vr) 2 

a, 



+ 



3$l(-2A fe /fc 2 ) + <5>\{-2a\ k /k 2 )} 



(34) 
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Here $>p(w) and are concerning with the integrals of the cutoff function, and defined by 



Tin 



1 r°° , „_iR i0) (x)-xRW(x 
dxx 



o 

l r°° 
o 



x 



+ R(°)(x) + 



w 



p ' 



— I <!xx n ~ l - Tv- 



Here, the cutoff function R(°\x) satisfies constraints given by Eq. (pi). In fJ], this is given by 



R(°\x) 



exp (x) — 1 

Now, we introduce the dimensionless Newton constant g k and cosmological constant X k , 



(35) 



g k = k 2 G k = k 2 Z N \G, X k = k 2 A 



If Eqs. (|33| ) and (^J) are expressed in terms of dimensionless couplings, we obtain 

P g = d t Qk = {1 + ri)g k , 

fix = d t \k = -(2 - rj)\k + g k B 3 (X k ,a). 

In this case, the anomalous dimension is given by 

g k Bi(X k ,a) 



n 



1 - 9kB2{X k ,a) 

Now, Bi(X k ,a), (t = 1,2,3) in Eqs. © and @ is defined by 

Bi(X k ,a) = ^-f3*l(-2A Jfc ) + 2*l(-2aA fc )-4*i(0) 



3vr 



B 2 {X k ,a) 



-15$|(-2A fc ) - 3(2a - l)*l(-2aA fc ) - 6*1(0) 
-— h$\(-2X k ) + 2$\(-2aX k ) 
-15*l(-2A fc ) - 3(2q - l)$ 2 (-2aX k 



B 3 (X k ,a) = i-[6^(-2Ajfe)+4^(-2aA fc )- 8*^(0) 
-7 ? {3$i(-2A fc ) + 2$i(-2aA fc )}" . 

These results are same as that of pf. If a = 1, these reproduce the results of 



(36) 
(37) 



(38) 



4 The UV NGFP of QG 



n 



4.1 Gauge dependence 

In below, we ignore the dimensionless cosmological constant. This approximation is reliable 
when the dimensionless cosmological constant is much smaller than the IR cutoff scale: <C k. 
This approximation is applicable if we are interested in the local structure of the Universe. 

On the FP, the scale invariance is preserved. Hence, if we denote g* as the FP of the 
dimensionless Newton constant, this satisfies 



Q = (2 + V *)g*. 



(39) 



From Eq. (]39|) , it is recognized that the candidates of the FP are g* = and rj* = —2. Here the 
former is the GFP and exists independently of a. The latter is the candidate of the NGFP. The 
condition if = —2 reads to 

9* = n 77T ■ ( 4 °) 



B 1 (a)-2B 2 (a) 



Here, Bi(a), (i = 1,2) is given by 



1 



fli(a) = — *}(0)-6(a + 3)*|(0) 

07T L 



B 2 (a) 



1 

6-7T 



5*}(0) -6(a + 2)*§(0) 



Now, to calculate $-(0) and &\(0) (i = 1,2), we use the cutoff function given by Eq. (|35|). 
Hence we have 

*i(o) = y, *2(o) = i, $l(o) = i, *|(o) = |. 



Therefore, Eq. © b ecomes 



* 27T / 

y 3 \ 



7T 2 -114 

54 



(41) 



As immediately recognize, Eq. ( |41| ) hcts the singularity at casing 

= (71- -114)/54. The existence 

of a singularity is not the problem, since this point is not the GFP. In the range a s i ng < a, g* 
has a positive value except for the limit a — > oo. On the other hand, in the range a < a s ing 5 9* 
has a negative value except for the limit a — > — oo. In the limit a — ► dboo, the NGFP merges 
to the GFP. This problem will be discussed in Sec. 5. The behavior of Eq. (41) is shown in 
Fig. 1 (a). In this figure, a vertical long dashed line corresponds to the singular point. A solid 
line is the positive NGFPs, and a dashed line is the negative NGFPs. In below, we consider 
only the positive NGFP. Hence, the shadowed region in this figure is ignored. 
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Fig. 1 : The gauge dependence of the UV NGFP of the dimensionless Newton constant in d = 4 (a), and the RG 
flows of it in a = 1 (b). 



The typical RG flows are shown in Fig. 1 (b), if a is fixed to unity. Under the present 
approximation given by Eq. (^), only the FPs are reliable. However, when a = 1, it is possible 
to study the behavior of the RG flows |], |(|. In this figure, the horizontal bold solid line, 
g* = 0.7152, corresponds to the UV NGFP (the circle in Fig. 1 (a)). This line separates the 
phase space into two regions; the strong coupling phase and the weak coupling phase. 

4.2 Cutoff function dependence 

As mentioned previously, constraints for the cutoff function are given by Eq. (^). Hence, any 
functions satisfying these conditions are applicable. Now, we slightly modify Eq. (^) as 

R? ) (x,s)= SX s>0. 
exp(sx) — 1 

Here, s parameterizes the profile of the cutoff function. As same as the previous subsection, the 
candidate of the UV NGFP is given by 

9* = — 7 7 7- ( 42 ) 

Bi(a, s) — 2B2{a, s) 

Here, Bi(a, s), (i = 1,2) is defined similarly to Bi(a), except for $-(0, s) and s), (i = 1,2) 
depending on s. 

If we numerically calculate the integral of ^(0, s) and &\(0, s) and insert these into Eq. (|4^), 
we get Fig. 2. In this figure, the surface represents the position of the UV NGFPs. The region 
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Fig. 2 : The gauge and cutoff function dependence of the UV NGFP of the dimensionless Newton constant. 

above the surface is the strong coupling phase, and that below the surface is the weak coupling 
phase. From this figure, it is confirmed that the UV NGFPs exist in a global range of a and s. 
The gauge dependence is same as the previous subsection. For the cutoff function dependence, it 
is recognized that the UV NGFPs merge to the GFPs in the limit s — > 0. The reason is expressed 
as follows. In the limit s —* 0, the integration of the cutoff function diverges, therefore, B{(a, 0) 
function diverges. Hence, the denominator in Eq. (^) goes to infinity, and g* goes to zero. 
However, there is no IR cutoff in the limit s — > 0. Hence, this limit is out of the applicability 
of the ERGE. If we take another type of the cutoff functions such as exp(— sx), same structure 
is observed. Therefore, the cutoff function dependence does not cause the disappearance of the 
UV NGFP and the change of the phase structure. 

5 Summary and discussion 

In this article, we consider the gauge and cutoff function dependence of the UV NGFP. For the 
cutoff function dependence, the UV NGFP strongly depends on the profile of the cutoff function. 
However, this dependence does not change the phase structure of pure QG. Though there are 
many cutoff functions, the phase structure will not be changed by this dependence. 

For the gauge dependence, the UV NGFP exists in a global range of a except for a = ±oo. 
This gauge, a = ±oo, is a bad gauge in this formulation. The disappearance of the UV NGFP 



14 



for this gauge is a serious problem. This is because, the phase structure of pure QG is changed 
in this gauge. This problem may be due to the treatment of the constant gauge parameter. 
Hence, if we will be able to treat the gauge parameter as the running gauge parameter, this 
problem may disappear. 

In this article, only the operators that are invariant under general coordinate transforma- 
tions are considered, because the gauge symmetry is preserved by the projection on = g^ u . 
However, if we improve this approximation to treat the running gauge parameter, the gauge 
symmetry is not manifestly maintained. Thus, the operators that can not preserve general 



coordinate transformations must be included in the functional space [12j. 

The other improvement is to formulate the ERGE without the gauge fixing. Recently, for 
pure SU(N) gauge theory, Tim R. Morris proposed the formulation of the gauge invariant ERGE 
|jl3|| . This formulation does not need the gauge fixing. Hence, if we can formulate the ERGE for 
QG in this formulation and show the existence of the UV NGFP, QG becomes an asymptotically 



safe theory and (non-perturbatively) renormalizable [14]. 

In this article, the Einstein-Hubert truncation is applied. However, for these two improve- 
ments, to accurately study the existence of the UV NGFP and the phase structure, we need the 
extension of the functional space. One possibility is to include the 72 2 -terms. The other is to use 
the method so called frame transformations |fL5fl . By this transformations, the higher derivative 
gravity is reduced to the Einstein gravity with the auxiliary tensor matter field. Hence, if we 
study this reduced theory, the phase structure of the higher derivative gravity will be clarified. 
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